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Abstract 

A quantum kinetic theory for correlated charged-particle systems in strong time-dependent electromagnetic 
fields is developed. Our approach is based on a systematic gauge-invariant nonequilibrium Green's functions 
formulation. Extending our previous analysis jlj we concentrate on the selfconsistent treatment of dynamical 
screening and electromagnetic fields which is applicable to arbitrary nonequilibrium situations. The resulting 
kinetic equation generalizes previous results to quantum plasmas with full dynamical screening and includes 
many-body effects. It is, in particular, applicable to the interaction of dense plasmas with strong electro- 
magnetic fields, including laser fields and x-rays. Furthermore, results for the modification of the plasma 
screening and the longitudinal field fluctuations due to the electromagnetic field are presented. 

1 Introduction 

With the progress in short-pulse laser technology H high intensity electromagnetic fields are 
becoming broadly available. In particular, they make it possible to create strongly correlated 
quantum plasmas under extreme nonequilibrium conditions which opens a broad range of 
applications, e.g. H. At the same time, optical techniques for time-resolved diagnostics 
are improving remarkably Q. These developments create the need for a quantum kinetic 
theory of dense nonideal plasmas in intense laser fields. 

Nonequilibrium properties of dense plasmas in which collisions are important are usually 
studied on the basis of kinetic equations of the Boltzmann type. However, in spite of their 
fundamental character, Boltzmann-like kinetic equations have a number of shortcomings, 
in particular in view of their application to dense plasmas in intense laser fields: 

i) they are valid only for times larger than the correlation (or collision) time r corr ~ , 

ii) they conserve only the mean kinetic energy instead of the sum of kinetic and potential 

energy, 

iii) they are valid only in the weak field limit since the corresponding collision integrals are 

independent of the electromagnetic field, 

iv) they are not applicable to high-frequency processes (fields), where uj > u p i, cf. i). 

Obviously, in the case of strong correlations, high-frequency electromagnetic fields and/or 
short-time phenomena generalizations are necessary. 

Generalized kinetic equations for correlated plasmas have been derived already in the 
60ies by Prigogine ||, Zwanzig ||, Kadanoff and Baym J?], H|, Balescu [H, 10 1, Silin [11] 



Klimontovich [12, [13, 14| and others. In recent years, the increasing interest in ultrafast 
processes has revived the theoretical activities, e.g. (15|, [R], [L7], |18|] , acompanied by progress 
in numerical solutions, e.g. [15, 19, p(I]; for textbook overviews, see [21, 22]. Furthermore, 
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kinetic equations for classical plasmas in high-frequency fields have been derived for the 
first time in papers of Silin (23|, f 



25]. Among other problems, he computed the high- 



frequency conductivity of a plasma. We mention that, for the weak-field limit, this problem 
has been studied by many authors, including Oberman et al. |2(3, and DuBois et al. 
An essential further development of the theory has been given by Klimontovich and 



co-workers |2 
phase space fl3 



J. Klimontovich used his powerful technique of second quantization in 
12| to investigate the density-density and microfield fluctuations in low and 
high-frequency fields. This allowed him to derive collision integrals for classical plasmas in 
strong fields which take into account dynamical screening and to derive a complete theory 
of transport processes [29, [HI]. A result of central importance is an expression for the 
collisional heating rate and the electron-ion collision frequency in strong fields in terms of 
the imaginary part of the inverse dielectric function line -1 [14]. Recently, expressions of the 
same form were derived again 32]. For a recent overview on the collision frequency in laser 



plasmas, we refer to Mulser et al. [33 



The above kinetic theories for plasmas in electromagnetic fields were limited to classical 
plasmas. A first extension to quantum plasmas was given by Silin and Uryupin More 
recently, a kinetic equation for dense quantum plasmas in strong static fields has been 
derived [[}(], [37j , whereas a systematic quantum kinetic theory for plasmas in strong fields of 
arbitrary time dependence was presented in ref. Q . There, electron-electron and electron- 
ion quantum collision integrals were derived within the static Born approximation (Landau 
collision integrals). 

In this paper, we extend this theory to the case of full dynamical screening for dense 
quantum plasmas under arbitrary nonequilibrium conditions by using the random phase 
approximation (polarization approximation) for the particle-particle scattering processes. 
Our approach is based on the nonequilibrium Green's functions formalism which allows for 
the most straightforward derivation and for an explicit solution of the gauge problem. We 
derive a kinetic equation which is a generalization of Klimontovich's classical result [p9| and, 
on the other hand, generalizes previous quantum results for the case of zero field [ 38 , [39|1 and 
static field [4C[. Furthermore, we derive results for the polarization and screening properties 
and longitudinal field fluctuations in a strong field of arbitrary time-dependence. 

2 Basic physical problems and definitions 

We consider the time evolution of a dense charged particle system under the influence of a 
strong time-dependent electromagnetic field and inter-particle correlations. 

2.1 Free particle motion 

It is instructive to recall first the motion of classical free charges in an external field. From 
integrating Newton's equation, m a dv a /dt = e a B{t), we obtain the velocity change of a 
particle with charge e a and mass m a in the field B(t) during a time interval [t',t], 



Av a (M') = — f cftE(t), 
m a Jt' 

and the field induced displacement 

Ar a (M') = — f di ! diBit) = — ! dt C diBit). 
m a Jf Jt' m a J t ' Ji 

In the equations above, we have dropped contributions from the acceleration and velocity 
at the initial moment t' since they are not related to the field. A further important quantity 
is the average kinetic energy 

fT dtvl(t), (3) 



(1) 



(2) 



2 T 
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where, in case of a periodic field, T is the oscillation period. Below we will also need the 
change of the relative velocity of a particle pair a, b gained in the field E(t) 



(4) 



Av ab (t,t') = Av«(t,0 - Av 6 (M') = (*L-*l) /*diE(i), 

\m a m b J J t , 

and the change of the inter-particle distance Ar ab 

Ar a6 (M') = Ar a (t,0 - Ar 6 (t,i') = ( — — ^ I dt I diB(t). (5) 

\m a m b J J f J i 

The most important special case is that of a harmonic time dependence 

E(t) = E cos ftt, (6) 

which leads to the following explicit results for the quantities introduced above 

Ava(M') = v° [sin ttt - sin M'] , (7) 
Ar a (t,t') = r° a [Q(t - t') sinQt + cos £lt — cos fit'] . (8) 

Here, we introduced two important quantities, the "quiver" velocity 

v° - ^% (9) 
m a \l 



and the so-called excursion amplitude 

o _ e a Eo 



m a £l 2 



(10) 



Furthermore, the cycle averaged (T = 2-ir/Q) kinetic energy gain of a charged particle (||) 
is the so-called ponderomotive energy 

2 p2 

pkin _ p pond ^q^O (i i \ 

" a 4m a ft 2 ' 1 j 



Obviously, the above results may be extended to electromagnetic fields with arbitrary 
time dependence, for example by expanding the field in terms of harmonic components. 
Nevertheless, it is useful to explicitly consider a second situation frequently encountered 
in modern applications: pulsed fields, e.g. those produced by femtosecond lasers. We will 
consider pulses of the following form 

E(t) = E p (t) cos Qt, E p (i) = 2E sin Q p t, < tt p t < it, (12) 

and E p = otherwise. Typically, S7 P <C S7, although modern femtosecond laser pulses may 
be as short as a few periods of the main frequency £1. For the field (12), we obtain the 
velocity change and displacement 

Ava(M') = -^vHsin^-sinttY], v°± = fi± = O p ± O, 

Ar a (t,0 = J^r° s [-n s (t-0cosn s t + sinn s t-sinnV] , 

with r° = e a Eo/m a f2 ±2 , whereas the change of relative velocity and two-particle distance, 
Av a j and Ar a 6 follow from Av a and Ar a by replacing by v^— and by r^— r^, 
respectively. 
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The above results trivially include the case of a time-independent electric field which is 
recovered by letting in Eq. (0) Q — > 0. The corresponding results are 

Av a (t, t') = ^E (t - 0, Ar a (t, t') = J— E (t - t') 2 - 

It is instructive to consider a number of parameters which characterize the state of the 
plasma, field strength, quantum properties etc: 

1. The field strength can be characterized by the ratio a = v^/v^a of the amplitude of 
the oscillation velocity (quiver velocity) v®, Eq. (0), to the thermal velocity vth,a = 
(fcT/m a )V2. 

2. The relative importance of the field and of particle-particle interaction is characterized 
by (5 = r®/rD, where is the amplitude of the field induced displacement (||), and 
rjj is the Debye radius (interaction range). 

3. The frequency of the field has to be compared to the eigenfrequencies of the plasma, 
most importrantly, the electron Langmuir (plasma) frequency, 7 = Q/u p i, which re- 
flects competition between field frequency and plasma density effects. 

4. The relevance of collisional processes depends on the ratio 5 = u/D,, where v is the 
total collision frequency of electrons in the plasma. 

5. The photon energy is characterized by its ratio to the thermal energy, hVt/kT. 

Modern lasers easily produce strong fields which satisfy the inequalities a ^> 1 and 3> 1. 
In high-frequency fields and/or plasmas of moderate density, 5 <C 1, which allows to treat 
collisions perturbatively, see below. 



2.2 Two— particle scattering 

Coulomb interaction between the charged carriers as well as quantum effects, obviously, may 
drastically modify the free particle behavior. Scattering of two particles with charges e a , e b 
in quantum states |ki) and |k2) on the Coulomb potential 

4Tre a e b h 2 

Vab{q) = ' ( 13 ) 

leads to a transfer of momentum q between them, so after the collision time t ~ t co \\ 
the particles are in momentum states |ki + q) and |k2 — q). While conventional kinetic 
approaches treat collision as instantaneous, t co n — > 0, this is not appropriate for correlated 
plasmas as well as in the presence of rapidly varying fields with f2 • t co \\ not being small. 
In this case, during the collision time, the scattering partners will be accelerated by the 
external field which is called intra-collisional field effect, which essentially modifies the 
scattering process. Using a quantum language, during the collision time, the particles may 
absorb photons of the electromagnetic field which is the familiar inverse bremsstrahlung, or 
re-emit them (bremsstrahlung). 



The kinetic treatment of two-particle scattering on the Coulomb potential (13) leads 
to the well-known divergencies at short and long wavelengths. While the first is naturally 
cured by a quantum theoretical approach, the origin of the latter is the long range of the 
Coulomb interaction. The familiar solution lies in the replacement of the bare Coulomb 
potential ( |l3| ) by a screened one 

V ab (q)^V: b (co, q ) = ^^- y (14) 
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where is the retarded dielectric function. This is not only of fundamental interest but has 
also practical relevance. The dielectric function includes collective plasma oscillations and 
instabilities which, especially in nonequilibrium situations, may strongly enhance scattering, 
transport and energy exchange with the electromagnetic field (anomalous transport). 

The simplest approximation for the dielectric function is the random phase approxima- 
tion (RPA) being the quantum generalization of the Vlasov dielectric function 

e R (u,q;t) = 1 - £ V aa (q) U R (u, q; t), (15) 

a 

U R (ua-t) - 1 f d3k /q(k;t)-/ a (k + q;t) 

where e a denotes the single-particle energy. Interestingly, this result was derived by Klimon- 
tovich and Silin [41, [4^] two years before Lindhard p^ ]. While this approximation is appli- 
cable to nonequilibrium situations in which the Wigner distributions f a are weakly time- 
dependent (when the time scale t is much longer than 2ir/u>), on short times or for fast 
processes, such as in high frequency fields, generalizations are necessary. Such generaliza- 
tions avoid the assumption of separation of the two time scales and lead to an explicit 
dependence of the functions II, e and V s on two times. Moreover, the presence of an electro- 
magnetic field may be expected to modify the dielectric and screening properties, leading to 
a very complex problem of coupled particle, screening and field dynamics. The appropriate 
theoretical concept to tackle this problem is provided by quantum field theory. 



3 Quantum field theoretical approach to the dynamics of plasmas in 
electromagnetic fields 

Numerous concepts have been developed to describe the mentioned above dynamics of parti- 
cles and fields. Among them, the most systematic and powerful is the theory of nonequilib- 
rium Green's functions. It is based on the method of relativistic quantum field theory, where 
charged particles and the longitudinal and transverse electromagnetic field are described on 



equal footing by field operators [44, 45 1. From the equations of motion for the field op- 
erators - the Dirac equation and Maxwell's equations, one can derive equations of motion 
for all quantities of interest. Among them, the most important are two-time correlation 
functions (Green's functions) which allow for systematic and far-reaching generalizations 
of traditional kinetic theory. In this paper, we focus on nonrelativistic particle dynamics 
and start our derivations from the familiar equations of motion for the particle correlation 
functions g > and g < while the electromagnetic field is treated classically. 



3.1 KadanofF Baym Equations 

The field theoretical description of plasmas is based on the creation and annihilation oper- 
ators ift anc [ ip [gg| which are defined to guarantee the spin statistics theorem, 

Va(l)^(2) T^(2)Va(l) = ^(1)^(2)T^(2)^(1)=0, 
^(1)^(2)^(2)^(1) = $(1-2) ^6, 

where t\ = ti has been assumed. The upper (lower) sign refers to bosons (fermions), 
1 = (ri,ti,sf), and a labels the particle species. Below, we will drop the spin index and 
assume fermions. The nonequilibrium state of a correlated plasma is described by the two- 
time correlation functions which are statistical averages (with the initial density operator 
of the system) of field operator products 

g>(l, 1') = 1(^(1)^(1')} , <?<(!, 1') = -4(^(1')^(1)), (17) 
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where g> and g< are, in nonequilibrium, independent from one another. They contain 
the complete dynamical and statistical information. The latter follows from their elements 
along the time diagonal: the one-particle density matrix is immediately obtained from the 
function g K according to 



/ a (r 1 ,r / l! t) = -^<(l,l')| tl=ti) 



(18) 



whereas the dynamical information (e.g. the single-particle spectrum and the correlations) 
follows from the function values across the diagonal in the t\ — t^-plane, in particular, from 
the spectral function o(l, 1'), 



a(l, 1') ee ih{g>(l, 1') - g<(l, l')} = ih{g*(l, l') - g A (l, l')} , 



(19) 



where g R l A are the retarded and advanced Green's functions, defined below in Eq. (25). 
In the following, it will often be convenient to use microscopic and macroscopic time and 
space variables being defined as 



r = n - r 1: 

T = t 1 -t' 1 , 



R = (n + r[)/2, 
t=(h + t[)/2 . 



(20) 
(21) 



In particular, in cases where the microscopic variables vary on much smaller scales than the 
macroscopic ones, it is advantageous to perform a Fourier transformation with respect to r 
and/or r which leads to the frequency and momentum variables u and p, respectively. In 
particular, Eq. (|T8|) then yields the familiar Wigner distribution function 



/ a (p,R,t) = -z'%<(p,R;ti,ti)| ti=< v = 



(22) 



The time evolution of the correlation functions in an electromagnetic field is determined 
by the Kadanoff-Baym equations 0, 47 1 



■ih 



d l_ 

dti 2m a 

ti 



Vi-^A(l)) -^(1) g%(l,l')- Jdri^il^ti^frtx,!') 
= Tdl [S>(1,I)-E<(1,I)] ^(1,10- f 1 dLVf(l,I) [g>(l,l')-g<{l,l')], (23) 

J t J t 

which have to be fulfilled together with the adjoint equations. Here, to denotes the initial 
time where the system is assumed to be uncorrelated (otherwise, the equations have to be 
supplemented with an initial correlation contribution to S a , cf. J48|] ). is the Hartree- 
Fock selfenergy (mean-field energy with exchange), 



^HF, 



-ih5{t\ 



t V ) Y,{f d ^ V ab{^l ~ r2)^(22 + ) - SafiV^Tx - ^{ll') 



) , (24) 



and Sa are the correlation selfenergies which will be discussed below. 

For the following derivations, it is useful to introduce, in addition, the retarded and 
advanced Green's functions 



gV A (l, 1') = ±6[±(f! - t[)] {g>(l, 1') - g<(l, l')} , 
which obey the simpler equations 

\f /A (M') 



(25) 



., d 1 h e a 
ih- -Vi A(l 

9*1 2m a \ i c y ■ 

- J d2SfM(l,2)^M(2,l') = 6(1-1'). 



(26) 
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In equations (|23|) and (^6|), the electromagnetic field is given by the vector and scalar 
potentials A and 4> an d will be treated classically. A is the full vector potential (external 
plus induced) which obeys Maxwell's equations, whereas <p is understood as to be due to 
external sources only, the induced longitudinal field is fully accounted for in the screened 
Coulomb potential V s which enters the selfenergies S^ - and , see below. 

Although one can directly analyze and solve the two-time Kadanoff-Baym equations 



(|23|), e.g. [19|, for an overview, see |22|| , it is easier to consider the kinetic equation for the 
Wigner distribution function (122|), which we will be concerned with here. This equation is 
immediately obtained from the equal time limit, t\ = t[ = t, of Eq. (|23|) plus its adjoint. 
Introducing further the variables R and r, Eqs. (|20| ) and (|2l|), we obtain, after Fourier 
transformation with respect to r, for the spatially homogeneous case, 

^/a(p, t) = -2 Re jf di{^> (p; t, t) g< (p; i, t) - £< (p; t, i) g> (p; t, t) } . (27) 

This is an exact equation and, therefore, well suited for deriving generalized kinetic equa- 
tions. However, this equation is not closed yet since it contains under the collision integral 
functions depending on two times. Therefore, to obtain explicit expressions for the collision 
integral, one has: 

1. to find appropriate approximations for the self energy. For this, the Green's functions 
approach provides powerful approximation schemes based on Feynman diagrams which 
allow for a very systematic development of the theory. Here, we are interested in the 
plasma dynamics with screening effects properly included. Therefore, the appropriate 
choice for the selfenergy will be the random phase approximation (RPA); 

2. to express the correlation functions as functionals of the Wigner functions f a (re- 
construction problem). This problem can be solved approximately on the ba sis o f the 
generalized Kadanoff-Baym ansatz (GKBA) of Lipavsky et al. fl4S[| , see Sec. |3.3| . 

3.2 Gauge invariant Green's functions 

It is well known that the electromagnetic field can be introduced in various ways (gauges) 
which may lead to essentially different explicit forms of the resulting kinetic equations. 
Although alternative derivations are successfully applied too, gauge invariance becomes a 
particular problem if the resulting kinetic equations are treated by means of approximations, 
such as retardation or gradient expansions. A critical issue is that the result of these 



approximations maybe essentially different in different gauges, see e.g. (21] for examples. 
To avoid these difficulties, we will formulate the theory in terms of correlation functions 
which are made explicitly gauge-invariant. While the main results have been presented in 
Ref. 0], here we provide some additional details. 

In this section, we use a co-variant 4-vector notation as it makes the following transfor- 
mations more compact and symmetric. The corresponding definitions are 

A^ = (c(f>, A), x il = (ct, r), X^ = (ct, R), 

and the convention = ao&o — ab is being used. 

One readily proofs that the Kadanoff-Baym equations ( |23| ) remain covariant under gauge 
transformations, i.e., under the following transformations of the potentials and field opera- 
tors 

A'^x) = A^x) - d^x(x) <(*) = etex^Mx), (28) 
The corresponding gauge transform of the Green's functions leads to 

a' a (x,X)=er-fH^)-x(^)] ga (x,X). 
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Following an idea of Fujita []50|], we now introduce a gauge-invariant Green's function g(k, X) 
which is given by the modified Fourier transform 



X 



(27T) 



4 exp < i 



dX x. 



(29) 



where use has been made of the identity 



X(X + 



x{x 



x 



2 d\ 4rX (X + Xx) = I ' dX \ ( X + Xx i . 

I dX 

2 



Indeed, one readily confirms that under any gauge transform (28), the phase factors cancel, 
and</(k,X)=5(k,X), P. 

In the following, we focus on spatially homogeneous electric fields and use the vector 
potential gauge 



A) = 4> = 0; A 
In this case, relation (p9|) simplifies to 



dtE(t). 



(30) 



<?o(k, u; R, t) = / (frdrexp 



i lot — — r 



V 



5a (r,r;R,t), (31) 



what means that the gauge-invariant Green's function g(k) follows from the Wigner trans- 
formed function g a (p) by replacing the canonical momentum p by the gauge-invariant 
kinematic momentum k according to 



ki a 
H 

c 



dt 



, A(t') 



(32) 



*-5 



Let us illustrate this for the examples studied in Sec. For a harmonic electric field given 
by Eq. (||), the vector potential and the momentum relation become, according to Eq. (|30|), 

Am cE ° • m i , 2E ■ n+ ■ fir 

Am = — — smstt, p = kH — -r smiltsin — . 

\l T 14 2 



(33) 



Similarly, for a pulsed field of the form Eq. (|i"2|), the result for the vector potential and the 
momenta is 



A(t) = cE ^ 



cos Sl s t 



=± 



2Eo \ -< 1 ... 
p = k H > — - cos i i i sin 

1 T ^ OS 



Finally, for a static field we obtain, 

A(i) = -cE Q i, 



p = k + e a E t. 



For the derivations below, we will need the gauge invariant Fourier transform of the 
convolution of two functions which, in the homogeneous case, is given by 



/(ri-ri; *!,*!)= / ^fS(n-f; fc, t) C(f - r' i; t, t^). 



(34) 
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After straightforward manipulations which involve the back transform of (0), we arrive at 
I(k;ti,fi) = 



1 dtB 









h , „ A(t") e a f h , „ A(t") 



*i - A 



tx-r 



x C 



, e a f* 1 , „ A(t") e a /■* , A(t") - , 
k + — / dt" T -±-L - -± I dt" t, t[ 



c J f > 



h - 1[ c y t / 



(35) 



In particular, the derivation of the collision integral in the kinetic equation for the Wigner 
function, requires the equal-time limit of this expression, t\ = t± = t, 



J(k; t) 



dtB 

x C 



k+^A(t)-— / dt 



t - 1 



tj 



c w c k t-t 



■ t,t 



(36) 



where B and C will be replaced by g^ and £<. Notice that in this case the momentum 
arguments of B and C are equal. To simplify the notation below, we introduce the field 
induced momentum shift 



which has the important property 



VA( t )-A( 



t-t' 



(37) 



K A (t,t') -Kf(t',t) = ^{A(i)- A(t')} = -ea f dt"V{t") = Q a (t,t'), (38) 

where Q a is nothing but minus the momentum gain of a free particle in the field, Q a (i, i') = 
-m a Av a (t,i'), cf. Eq. (jl]). Another important relation follows from multiplication by the 
time interval: 

— K A (t,t')-(t-t') = -R a (t,t'), (39) 
m a 

where R a is just the field induced displacement of a free particle, Eq. i.e. R a (i,i') = 
Ar a (M'). 

Definition (|37|) allows us to rewrite Eq. (|3"6|) as 



I(k;t) = ydtB [k + K^(t,t); t,t] C[k + K^(t,t); t,t] . (40) 
3.3 Gauge invariant propagator. Generalized Kadanoff Baym ansatz 



As noted in Sec. ^lj, for the derivation of the collision integral in the kinetic equation, we 
need to express the functions g > and g < in terms of the Wigner function. In addition, such 
a reconstruction ansatz involves the retarded and advanced Green's functions g R / A (t,t') for 
which suitable expressions have to be found. We first determine these quantities for free 
particles in an electromagnetic field which allows to simplify Eq. (26) to 



gV A (p;t ll t' 1 ) = 5(t 1 -t' 1 ) 
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which is solved immediately by 



9a(p;r,t) = --9(r)exp 



.» / dt f [p-^A(t')] 2 /2m a 



(41) 



and is obtained from the symmetry relation g A (p;r,t) = [g^(p; — r, t)]*. From this 
result, we can calculate the spectral function a(t,t ; ), Eq. (|i~9|), 



a a (p;T,t) = exp 



.« / rff [ p _ ^A(t')] 2 /2m a 



(42) 



Obviously, the results (f|l]) and (|42|) are gauge-dependent since ga^ A and a a are functions 
of the canonic momentum p. But one can easily obtain the corresponding gauge-invariant 
results by applying the transform (|29|), with the result 



where 5 a (A;r, t) 



.Ie(r)e-*L^ T+5o(A;T ' t) 
h 



2m a c 2 



2 1 / /•*+§ 



dt'A 2 (t') + - 



dt'A(t') 



(43) 
(44) 



This result has a simple physical interpretation. For a free particle without field, the spectral 
function shows free undamped oscillations along r (i.e. perpendicular to the time diagonal) 
with the one-particle energy e a (/c) — k /2ni a , and its Fourier transform is 



ai ree (k;u,,t) 



5[hjj - e a (k)]. 



(45) 



This clearly underlines the meaning of the functions g R ^ A and, in particular, the spectral 
function a - they contain the full information on the single-particle energy spectrum. Fur- 
thermore, in a correlated system, the single-particle spectrum is affected by interactions 
with other particles. This leads to a shift of the oscillation frequency and to damping of 
the oscillations, i.e. to finite life time effects, and it is reasonable to call the corresponding 
single-particle excitations quasi-particles. On the other hand, the result (43) reflects the 
influence of an electromagnetic field on the particle spectrum, while correlation effects have 
been neglected. Eq. (|4^) shows that the field causes a time-dependent shift of the single- 
particle energy which, obviously, reflects the well-known fact that the proper eigenstates 
of the system contain the electromagnetic field and are given by Volkov states [51]. The 
spectrum may even contain additional peaks which becomes particularly transparent in the 
limiting case of a harmonic time dependence: For the field (p), the time integrations in S 
can be performed, and simple trigonometric relations lead to [ j2[ 



5 a (A;r,t) 



eP° nd r 



sin VLt cos 2Vtt 

1 nr + 



1 sin 2 Sit sin 2 ^ 



(46) 



where Ea° nd is the ponderomotive potential which was introduced in Eq. ([Tl]). The first 
term in the brackets leads to a shift of the single-particle energy, the average kinetic energy 
of the particles increased by Ea° nd . The remaining terms modify the s pec trum qualitatively 
giving rise to additional peaks which are related to photon sidebands [52|. 
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Now we turn to the solution of the reconstruction problem. The simplest solution is the 
common Kadanoff-Baym ansatz, 

±ihg$(p;«;,t) = 4 ree (p; u,, t)/^(p; i), (47) 

where / < = / and / > = 1 — /, and the upper (lower) sign refers to g > (g < ). Indeed, 
the two-time functions g^^t\A\) and are now expressed in terms of one-time 

Wigner distribution functions and a known spectral function. However, due to the expected 
retardation effects, this ansatz is not applicable here. As mentioned above, a more general 
solution which properly takes into account retardation (memory) effects is the generalized 
Kadanoff-Baym ansatz proposed by Lipavsky et al. p9[ which reads 



gfiv-hA) = iKg?(p;hA)9%(p-AA) -*^(p;*i.*i)^(p;*i»*i). ( 48 ) 



where for the functions on the time diagonal ±iHg£(p; t\, tf) = /<(p; ti), cf. Eq. fl22|). 
Within the quasiparticle approximation and with static selfenergies, Eq. (^) is exact. In 
more complex situations, it is an approximation to the exact reconstruction solution, which 
has prooved extremely successful in many applications. In particular, it has been used 

for more general selfenergies and also with more general propagators g a , e.g. [21|. We, 
therefore will use this ansatz below. 

Eq. (fig 1 ) is written in terms of the momentum p and is, therefore, gauge-dependent. To 
transform this relation into a gauge-invariant form, we use its coordinate representation, 

±g*t(r 1 -r' 1 ;t 1 ,t' 1 ) = j dr g* (n - v-txA) f£(? ~ 

- J dvf^-r^gair-r'iMA), (49) 



and apply the gauge-invariant Fourier transform ( pip together with the back transforms of 
ga^ A and f£~ which leads to the gauge-invariant generalization of the GKBA 

±g^MA)=9a^MA)f^ [ k - K a -ff- L k-K^(M');*i] ^(Mi,<i), (50) 

where the definition (^) for has been used. As in the field-free case, the first term is 
nonzero only for t\ > t 1 and the second in the opposite case. Notice the difference of the 
time arguments in the two distribution functions. 

4 General kinetic equation for quantum particles including screening 
and electromagnetic fields 



We now come back to the time-diagonal limit of the Kadanoff-Baym equations, cf. Eqs. (p 
and (p7[), and derive the quantum kinetic equation for a plasma in a laser field thereby fully 
taking into account dynamical screening. Again, it is advantageous to derive this equation 
for the gauge-invariant Wigner distribution. To this end, we take the Fourier transform 
(31) of the time-diagonal Kadanoff-Baym equation 

^- t f a (K,t) +e a E(t) ■ VMKA = "2Re^ dt{^>g< - S< 5 >} = I a (k a ,t), (51) 

where the full arguments are, according to the convolution relation ([5]), given by 

Xfgf = S| [k a + Kf(t, t); t, i] gf [k a + K£(t, t);i, t] . 

This expression is valid for arbitrary approximations for the selfenergies S > and S < . In our 
previous paper Q, we used the simple static Born approximation. Here, we are interested 
in a fully selfconsistent inclusion of dynamical screening, so the appropriate choice is the 
random phase approximation (RPA). 
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4.1 Random phase approximation 

Starting from the familiar expression in coordinate representation, application of the gauge- 
invariant Fourier transform (|3l]) straightforwardly leads to the following gauge-invariant 
result 

£|(k;M')=^/ sl(k-q;M')^(q; tj), (52) 
which transforms the collision integral of Eq,([5l"l) into 

I a (k a , t) = -2 Re J jf di | 5 > [k a - q + K A (t, t);t, i] 

x ihV^fattfgZ [k + K A (M~);t,i] - [>^-><]j (53) 

In Eq. (]52] ) we introduced the correlation functions of the screened potential (plasmon 
correlation functions) V>,V<, which contain the whole screening problem and are directly 
related to the correlation function of the longitudinal field fluctuations (microfield fluctua- 
tions) m 

4f WWW))* = \ [i^>(q; t, t') + i/rt£<(q; t, t')] . (54) 
q z 

While in the classical case, the contributions from V > and V < are equal, in the quantum 
case a symmetrization is useful which is indicated by the bar over the fluctuation term. V s > 
and V s < can be related to the retarded and advanced screened potentials via the optical 
theorem 



v^( q; tl M) = Y. / 1 dt * I * ^(q; h,ta) nl(q; t 3 , t 4 )v; s b A (q; u, t 2 ), (55) 

„ J to J to 



dt 3 
1 to J t 

where and V A obey the following equation of motion (Dyson equation) 



V: b R M (q; t, t') = V a MW ~ t') + V Kc(q) /" ^ n£/ A (q; t, i)V c f M (q; t, t'), (56) 

and y sR / A are related to the nonequilibrium inverse dielectric function according to 

Vj* M (q; f , t') = VM [e fl / A (q; t, t')] . (57) 

In the above equations, V^>(q) is the bare Coulomb potential (jl^) and U R / A the retarded and 
advanced longitudinal polarization functions (plasmon selfenergies) . To close this system of 
equations, the polarization functions have to be expressed in terms of the particle correlation 
functions for which the simplest approximation is provided by the RPA, 

n|(q;ti,t 2 ) = -ih J -^^ E gf(k b + q;t 1 ,t 2 )gf(k b ;t2,t 1 ), (58) 
<(q;M') = e(t-t'){n>(qiM')-^(q;M')}, (59) 

and H A follows from the relation lTd(q; i, t') = [IlS(q; t' , t)]*. This set of equations com- 
pletely defines the non-Markovian polarization approximation (RPA) for a quantum plasma 
in a strong transverse field. 
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4.2 Application of the gauge invariant GKBA 

What is left to obtain a closed expression for the collision integral in equation (51) is to 
apply the gauge-invariant GKBA (50) together with the free-particle approximation (p 
to all two-time functions. This leads to the following results for the optical theorem: 



Vi*(q;ti,fe) 



d?k Z"* 1 f t2 
j^j^ jf dt 3 ^ dUV* b R (<i; t u t 3 ) V b s a A (q; i 4) t 2 ) 



e k+q 



-el){t 3 -U)-qR b (t 3 ,t 4 ) 



9(t 3 -U)e~^' 
xff [k + q + Q b (t 3 ,U);t A ] /f [k + Q 6 (t 3 , t 4 ); t 4 ] 

+e(t 4 - 1 3 ) e -M(4 +q - e k)(*3-*4)+qR i ,(t4,t 3 )] 

x/^ [k + q + Q 6 (t 4 , t 3 ); t 3 ] /f [k + Q 6 (i 4 , t 3 ); t 3 ] 



(60) 



Similarly, the result for U R can be transformed to 



n y (27rft) d 

{/a [k + Q a (t, t'); t 1 ) - f a [k + q + Q a (t, t');t'] } , 



(61) 



where the momentum shift Q a and field-induced displacement R a were defined above in 
Eqs. ( |38l) and Eqs. (|39|), respectively. In the absence of the electromagnetic field, (Q a — > 
0,R a — > 0), Eq. (|6l] ) reduces to the well-known nonequilibrium RPA-polarization function. 
The effect of the field is two-fold: first, it introduces an additional retardation Q a in the 
distributions (intra-collisional field effect) and second, it leads to a modification of the 
one-particle energies in the exponent given by R a which we discussed in detail in ref. fi]]. 

We now can transform the collision integral Eq. ([53]) by applying the GKBA (^) to 
g > ,g < , and using the result for V s> and V s< , Eq. (|60|). After straightforward calculations, 
we obtain the final result 



Ia(K,t) 



dk b dk a dk b 
1 1 2 ^ ^ J (2vr^) 6 



5(k a + k b - k a - k, 



,) fdt 
J t 



e -i[(6 a - ea )(t-t)+(k a -k a )R (t,t)] f> [k a + Q a ( t ,t);t] f< [K + Q a (t,t);t\ 



Jto Jto 



duv sR (K-K;h,t 3 )V SA (K-K;U,t 2 ) 



sA, 



9(t 3 - £ 4 ) e -l[( e 6- e 6)(*3-*4)+(k i ,-k i ,)R i) (t3,t4)] 

Xf> [k b + Qb(t3,U);t 4 ] / 6 < [k6 + Q6(t3,t4);t4] 
+9(i 4 _ t 3 ) e -|[(e6-e6)(t3-t4)-(k 6 -ki ) )R i) (t4,t3)] 



<] .(62) 



This is the general non-Markovian collision integral for a homogeneous weakly coupled dy- 
namically screened plasma in an electromagnetic field. It is a generalization of numerous 
results which were previously obtained by various authors, including our static screening 
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result for strong time-dependent field s pj and the result of Silin and Uryupin [ 35 1 , the RPA 
result for a static field of Morawetz [ |40| ] and the field- free RPA results of Kuznetsov [38] 
and Haug and Ell Furthermore, it generalizes previous results obtained for classical 

plasmas by Silin, Oberman et al., Klimontovich and others p3| , [27], ff3j. In particular, as we 
will see below, the classical dynamical screening result of Klimontovich and Puchkov 2£] is 
straightforwardly recovered from the collision integral ( p2] ) . This collision integral is the ba- 
sis for computing electron-ion collision frequencies, plasma heating inverse bremsstrahlung 
effects etc., thereby fully taking into account dynamical screening, plasma instabilities and 
anomalous transport 



Despite the complicated structure of the integral (62), a direct numerical integration of 
the kinetic equation ([5l]) appears to be within reach, as recently solutions of non-Markovian 
RPA-type equations (without longitudinal field) for semiconductors have been reported 
On the other hand, to gain deeper insight in the physical processes contained 



[20 



in the collision term (^), it is useful to consider analytical simplifications. 



4.3 High frequency fields. Silin ansatz 

If the collision frequency is low compared to the oscillation frequency of the field, i.e. if the 
parameter 5 <C 1, one may follow an idea of Silin |2^j and solve the kinetic equation ( |5l|) 
with a perturbation ansatz f a = fa+ fa where f% obeys the collisionless equation 

d 



-/ a °(k,t) + e a E(t).V k / a °(k,t)=0 



with the solution /„(k, t) = f, 



of) 



(63) 



where f a o is an arbitrary function depending on the initial conditions. The equation for f\ 
reads 



_/i(k,t) + e a E(t) • V k /i(k,i) = I a (k,t) 
where J^k, *)=!„{/«-► f a0 [k + ^A(i)] } . 



(64) 



With this scheme, there follow essential simplifications of the above results because it effec- 
tively eliminates the time retardation of the distribution functions in the collision integrals 
|p6fl . Indeed, one easily verifies that the arguments of the distributions which appear in the 
formulas above, now become 



f a [k + Q a (t,t');t'] — >/, 



«o 



k+^A(t) 

c 



(65) 



and do not depend on the time t' anymore. This simplification allows to compute the 
transport, screening and fluctuation properties quite efficiently. 

We first consider the modification of the longitudinal polarization IL R , Straightforward 
transformations of Eq. (^) including a change of the momentum integration variable lead 
to 



n£(q;t,0 

whith nf (q;r) 



■ / dtA(t) 



O (q;*-0 



(66) 



--6(r) 



d 3 k 



e-n^-^r {/ (k)-/ (k + q)}. (67) 



Using this result and the adiabatic approximation (neglecting the ion contribution to the 
polarization), we obtain from the Dyson equation, Eq. (|56[), for the retarded screened po- 
tential 



^f(q;M') 



m a ch 



J dtA(t) 



Kf,o(q;*-0 



(68) 
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where again, Vj£ Q denotes the screened potential in the zero field limit and without retarda- 
tion in the distribution functions. Simplifications are possible also for the collision integrals. 
In particular, we obtain for the electron-ion scattering term, 



c 

' "I 

CO 



dn I dr 2 — _ - e ^ {/ e0 (k - q) - / e0 (k)} , (69) 



r° , 1 1 /_ diA(f) 

'e^(q,Ti) e^(q,r 2 ) 

where l/e^ A are the field-free inverse dielectric functions. Eq. (|69| ) is the generalization of 
Klimontovich's result [14| to quantum plasmas and. 

Finally, we consider the simplifications introduced by the ansatz ( |64|) to the field fluctu- 
ations. The quantity of central importance is the correlation function of the temporal and 
spatial microfield fluctuations (5E5E), Eq. (|||). A lengthy but straightforward calculation 
leads to the following result 

(4vr) 2 ^/W") f d 3 k 



UttY l ^ J MA W f 

<aE(ti)fl3(<s,)) q = *» J 



(2-kK) 



, , , 1^ -^) 

e e 1 7 e 7^2 iM* + 0)[1 " /d)(k)] + / e0 (k)[l - / e0 (k + q)]} 



g ft L c k+q c k 

c e — 

e k+q c k 



+ e 2 e -iK +q -L](*i-*2) / dn 



* I J- /?/_ 

tX-TJ 



J — (X 



dr 2 ^ _ N e e w *2-^ [/ i0 (k + q) + / i0 (k)] 



where the first term in parantheses (second line) is the electron contribution, and the second 
(third and fourth line) results from the ions. Again, this is a generalization of Klimontovich's 
remarkable result [14] who considered the classical limit and the equal time fluctuations, 
t\ = t2- Our result fully includes the two-time fluctuations which are directly measurable 
quantities. From the above fluctuation spectrum, all major observables of dense quantum 
plasmas in a strong laser field can be computed. The corresponding analysis will be presented 
in a forthcoming paper. 



5 Discussion 



In this paper, we have presented a gauge-invariant derivation of the quantum kinetic equa- 
tion for dense plasmas in a laser field. Our main result, Eq. (|6^), generalizes previous 
work to quantum systems. This equation can be used to calculate the transport properties 
of a dense plasma in a laser field on arbitrary time scales, i.e., over the whole frequency 
range. The use of the random phase approximation allows for a highly consistent treatment 
of the combined effect of internal longitudinal fields (dynamical screening) and transverse 
electromagnetic fields, including intense laser pulses. In particular, it allows to investigate 
the influence of the electromagnetic field on the two-particle scattering process and the 
screening properties of the plasma and on the screening buildup in the presence of a strong 
field. 

Besides, the presented gauge-invariant approach is completely general and can be ex- 
tended straightforwardly to more complex situations, including strong coupling effects, 
bound states, impact and field ionization. Moreover, it can be directly generalized to rela- 
tivistic systems und ultra-intense fields. 
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